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The origin of the neutrino mass, dark matter (DM), and dark energy (DE) are among the most 
challenging problems of fundamental physics. We address these questions from analyses of the 
models where the neutrino mass is generated via Yukawa coupling to the quintessence field which 
represents the DE. It has been shown in a recent work on the toy model with a single Dirac 
fermion coupled to the quintessence that by choosing parameters of the DE potential to match 
the present DE density, the model allows to lock the neutrino mass at m ~ 10~ 2 eV and yields 
consistent estimates for other parameters of the Universe. To include the DM component into 
this framework we propose here to add the right-handed Majorana term(s) into Lagrangian with 
the quintessence-generated mass. As a results, the DE field is responsible for generation of the 
masses of the light active Majorana neutrinos as well as of the heavy sterile neutrinos. The latter 
are natural DM candidates. 
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1. Introduction 

It appears now that consensus has being reached on our understanding of the composition 
of the Universe and fundamental issues demanding extensions of the standard model of particle 
physics. Roughly 75% of the Universe is comprised of dark energy (DE), about 22% of it is 
comprised of dark matter (DM), and only approximately 3% corresponds to the contribution of 
the conventional (baryonic) matter. The present stage evolution of the Universe is governed by the 
dominant DE contribution, and the Universe experiences an accelerating expansion. The nature of 
the both DE and DM is still unknown, and this is probably one of the greatest mysteries of modern 
physics [JXJ] _ The origin of neutrino mass is also a fundamental problem for particle physics as well 
as for cosmology [Q, Si. Current upper limits on the sum of neutrino masses from cosmological 
observations are of the order of leV, while neutrino oscillations give a lower bound of roughly 
0.01 eV [gp. Even if at least ad hoc, minimalistically and intellectually unsatisfactory, the DE 
can be accounted for by adding a cosmological term to fit current observations, the existence of 
neutrino masses and DM particles definitely requires revisions of the standard model. There is a 
vast literature and a great interest in these issues. Due to space constraints we address the reader to 
"a short list" [|l]-|6|] where more references, especially on earlier work could be found. 

In this paper we summarize and advance further our recent work \Uft. It stems from the earlier 
proposal by Fardon et al [§]. In the approach [Q] the origin of the neutrino mass is due to the dark 
energy in the form of quintessence. Both the DE density and the neutrino mass are determined by a 
single mass scale of the DE potential, and consistent estimates for other parameters of the Universe 
are obtained. 

Here we propose how to incorporate the "third element of the puzzle", i.e, the DM component 
of the Universe within a common framework. For this end we analyze a quite straightforward 
extension of the "toy-model" of a single Dirac fermion considered in [0] by adding a right-handed 
Majorana fermion whose mass is due to the Yukawa interaction with the quintessence field. The 
Dirac-Majorana action for a single generation of neutrinos is known to be equivalent to two massive 
Majorana fermions. One of these, the (light) neutrino with mass m a decreasing in time can be 
identified with the conventional (active) neutrino, while the (heavy) sterile neutrino with growing 
mass m s is a natural candidate for the DM particle. The results and estimates we report below 
demonstrate that the theory is viable and consistent in unifying neutrino masses and the dark sector 
of the Universe. 



2. Toy Model: Dirac Fermion Coupled to the Ratra-Peebles Quintessence [|7]] 

The standard methods of general relativity and finite-temperature quantum field theory [^] 
were applied for fields in the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) Universe. The 
scale factor a(t) is governed by the Friedmann equations. The DE is modeled by the bosonic scalar 
field (quintessence) in the FLRW metric with the Euclidian action: 
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dx 



a(t) 3 d 3 x 



(2.1) 
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where (p = <p(x, t) and h = c = ks = 1. The analysis is restricted to the Ratra-Peebles quintessence 
potential: 

M a+4 

U(<p) = — , (2.2) 

with a > 0. In this toy model a single massless neutrino is mimicked by the the Dirac Euclidian 
action 

S F = J d% J a(t) 3 d 3 x y(x, t) (f-^ - l -J- v) y/(x, t). (2.3) 

The fermion and boson sectors of the Euclidian Lagrangian are interacting via the Yukawa coupling 
term: 

&yd = g(PYY ■ (2-4) 

Note that the Ratra-Peebles potential U((p) does not have a nontrivial minimum, so the generation 
of the fermion mass is due to the coupling between the quintessence and fermions resulting in 
a nonzero average of the bosonic field: m = g((p). (In the following the dimensionless Yukawa 
coupling is set g = 1.) Evaluation of the grand partition function in the saddle-point approximation 
[^] yields the density of the thermodynamic potential 

a((p) = u{(p)+a F {<p), (2.5) 

calculated at <p = (<p), where Q.p is the standard density of the thermodynamic potential of the free 
Dirac field [9]. The average of the scalar field is determined by minimization: 



dQ.((p) 



= U'((p) + gp ch = 0, (2.6) 



dip 

where we introduce the chiral fermionic condensate density: 

2m n tf-(Pm) 2 )l 
Pch = ToT / 7— i • (2-7) 

n 2 p 2 Jji m e z +\ 

Analysis of the extrema of the thermodynamic potential ( |2.6| ) reveals three phases of the coupled 
model. 



J ). Stable phase, T > T m : At high temperatures the equation (2.6) has two nontrivial solutions. 
The root indicated with a large dot in Fig. [T] (case a) corresponds to a global minimum of the 
potential. It is a thermodynamicaily stable state. In this phase CI < 0, so the pressure P = — Q. is 
positive. In the high-temperature phase the fermionic mass is small: 

m^M^V6ajY +2 oc T-ih , (2.8) 

The fermionic contribution is dominant and the thermodynamic potential behaves to leading or- 
der as that of the ultra-relativistic fermion gas. Two masses in this model demonstrate opposite 
temperature dependencies. The quintessence field with the conventionally defined mass 



2 d 2 U((p) 
m A — 



d(p 2 



(2-9) 

<P=(<P) 
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Figure 1: The dimensionless density of the thermodynamic potential Q.r = Q./M 4 as a function of the 
dimensionless parameter K = g<p/T. The thermodynamically stable solutions of Eq. ( |2.6[ ) indicated by the 
large dots correspond to the minima of the potential (curve a). The arrows indicate unstable solutions, 
corresponding to the maxima of the potential. The curve b is calculated at the temperature of metastability 
T m . At lower temperature T c the local minimum becomes an inflection point of the potential (curve d), and 
the system gets unstable. This is a point of spinodal decomposition. Adapted from 0. 



is "heavy" at high temperatures: 



m<sf 



a + 1 



T , M/T < 1 



(2.10) 



however its mass decreases together with the temperature. In contrary, the fermionic mass m 
monotonously increases with decreasing temperature. 

Another non-trivial root of (^^) corresponds to a thermodynamically unstable state (maximum 
of Q. indicated with an arrow in Fig. [[} There is a trivial third root (p = oo corresponding to the 
"doomsday" vacuum with a vanishing dark energy. In the high-temperature phase the doomsday 
vacuum CI = is a metastable state of the Universe. 

2) . Metastable phase, T m <T <T C : At the temperature of metastability T m the thermodynamic 
potential has two degenerate minima £l(((p)) = P(((p}) = H(oo) = 0. This is shown in Fig. [T] (case 
b). Below this temperature the two minima of the potential exchange their roles. The root indi- 
cated ny a large dot becomes a metastable state with Cl(((p)) > 0, i.e., with the negative pressure 
P(((p)) < 0, while the stable state of the system corresponds to the true stable vacuum of the Uni- 
verse £2(oo) = P(oo) = 0. See Fig. [j] (case c). The state in the local minimum Q.(((p)) is analogous to 
a metastable supercooled liquid. We disregard here the exponentially small sphaleron contributions 
of the fermions hopping from the metastable state £2((<p)) into the vacuum state £2(oo) = 0. 

3) . Phase transition at T = T c (spinodal decomposition): At the critical temperature 

M 

T c ^— (2.11) 
A c 

two extrema of the potential merge into an inflection point. This is shown in Fig. [T] (case d). Here 
we define 

v2 v _ v \ a+4 . 5 

V e 



3/2 



a + - 2 . 



(2-12) 
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From the viewpoint of equilibrium thermodynamics the model must undergo a first-order (discon- 
tinuous) phase transition and reach its thermodynamically stable (at T < T c ) phase corresponding to 
the doomsday vacuum 12(<p = oo) = ,p(<p = oo)=0. At the critical point the sound velocity and the 
compressibility vanish, thus the phase transition occurs through spinodal decomposition. During 
this transition the fermionic mass given at the critical point by 



and the quintessence mass 



m c K,^-M (2.13) 



nfyxs V / o(a+l)(^J M (2.14) 

both jump from the above values to their values in the vacuum state m = °° and = 0. 

It is quite a nontrivial problem to analyze the time evolution of the Universe towards the new 
equilibrium vacuum state. In our work [0] we abandoned completely the thermodynamic or kinetic 
considerations below T c and resort to the solution of the equation of motion 

to + 3Hto + ^= 0. (2.15) 

together with the Friedmann equations. The parameters found at the critical point were used as 
the initial conditions for the dynamics at T < T c . The numerical solution for the quintessence field 
(f){t) from the critical point to the present time oscillates quickly (with the period T ~ 10 Gyr) 
around the smooth ("mean value") solution (f>(t). The latter is found analytically as 

Relating the mean values with the physically relevant observable quantities, the key results are 
obtained analytically. To connect the only physical scale of the model M with the observables, we 
set the current density of the quintessence field coupled to the fermions to be equal to the observable 
value of the dark energy, i.e., roughly 3/4 of the critical density: 

Ppv.now = Pde ~ 7 • ^ ~ 31 • 1(T 12 eV 4 , (2.17) 
4 57TO 

and the combined baryonic and dark matter pM,now ~ | ■ ■ The exponent of the quintessence 
potential a is then the only parameter which can be varied. It is found that to get a realistic estimate 
for the present-time neutrino mass m[t) = (f>(t), one needs a small a. For instance, at a = 0.01 
m now « 0.27 eV. The relative energy densities Sl# = p#/p to t as functions of redshift (or temperature) 
are plotted in Fig. ||. There we also indicate the critical point parameterized by the redshift z c and 
the crossover point z* ■ The latter is defined as the redshift at which the Universe starts its late-time 
acceleration, i.e., where w to t = ^tot/Ptot = — |- For the present time we find w to t ~ — |. 

Being successful in relating the neutrino mass to the DE, this toy model have two major flaws. 
First, it does not give a clue about (heavy) DM particles. Secondly, a more serious problem is that 
in order to get a small neutrino mass we need a small (<p). This makes the mass too sensitive to the 
local temperature and thus the existence of the fundamentally identical particles in the Universe 
very problematic. As we will see in the next section the seesaw mechanism in the extended model 
allows to resolve these two problems on the same footing. 
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Figure 2: Relative energy densities plotted up to the current redshift (temperature, upper axis): Clyy - 
coupled DE and neutrino contribution; Q.y - radiation; Q.m - combined baryonic and dark matters. Parameter 
M = 2.39 • 10~ 3 eV (a = 0.01), chosen to fit the current densities, determines the critical point of the phase 
transition z c ~ 3.67. The crossover redshift z* ~ 0.83 corresponds to the point where the Universe starts its 
accelerating expansion. Adapted from P]. 

3. Extension of the toy model: Seesaw, Active and Sterile Majorana Neutrinos, 
Dark Energy 

To elucidate key points a theory needs to unify all three elements - neutrino, DM and DE 
within the same framework, we extend minimally the above toy model. For simplicity we take 



only one flavor of neutrinos. In addition to the Dirac mass term in the Lagrangian ( |2.4[ ) written in 
terms of the left- and right-handed components of the Dirac field y as Jzfyo = thj)Yr¥l + h.c, let 
us assume the existence of the right-handed Majorana neutrinos with the mass term m«: 

&m = ^m R %y R +h.c. , (3.1) 

where the charge conjugation is defined as = —iY 2 Yr- The extensions of the standard model 
to include the Majorana mass term have been known for quite a while (for reviews and original 



references, see, e.g., [|10|, [llj]), and more recent work pertinent to the present context in [|4j-|6j]). We 



take the standard assumption from the earlier work ([10, 11, f|-|6|] and more references there in), 
that the Dirac mass is due to the Yukawa coupling to the neutral Higgs scalar, i.e., mrj = go(H). 
We propose here that the origin of the Majorana neutrino mass is due to the coupling to the dark 
energy (quintessence), i.e., thr = g(<p). To the best of our knowledge this was not studied before. 
The proposal closest to the present one is put forward by Shaposhnikov and co-workers [|], ||] in 
the framework of the so-called vMSM extension of the standard model, where the right-handed 
neutrinos get their masses through the coupling to the inflaton field. 1 

Treatment of the fermionic action with the Dirac and Majorana terms is a well-known (seesaw) 
problem [[H], pp. The action can be brought to the diagonal form of the two Majorana fermions 
with positive masses: 

j(y/m 2 R +4m 2 D ±m R ^ . (3.2) 



m,, 



1 It is quite plausible that the inflaton and quintessence represent the same physical field analyzed at the different 
regimes of the Universe evolution. In this sense the current proposal might turn out to be equivalent to the one of |ph. 
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Taking the Dirac mass as a fixed parameter, we treat the model of two Majorana fermions coupled 
to the quintessence along the lines of the problem with a single Dirac fermion 0. In the saddle- 



point approximation we get the thermodynamic potential ( |2.5| ) with the fermionic contribution from 
two Majorana fields and minimization equation 

U'(q>) + - . — 1= (m s p ch (m, ) - m a p ch (m a ^j = , (3.3) 
' rrv\ + 4m?, 



Let us now set the physically acceptable range of the model's parameters. We set the quintessence 
Yukawa coupling g = 1 (this can always be done by rescaling of M [f7|]) and take the mass of the 
active light Majorana neutrino at present time m a = 10 eV. Since we don't know the Yukawa 
coupling go between the Dirac and Higgs fields ((H) = 174 GeV), we put nip in the range from the 
temperature of electroweak scale where the Higgs bosons acquire masses, i.e., mu ~ Tew ~ 10 2 
GeV to the scale of light leptons mp ~ 1 MeV. Then the current mass of the heavy (sterile) Majorana 
neutrino is actually given by the quintessence m s « niR « m D /m a and lies within 10 15 — 10 5 GeV. 



Matching then the current Ratra-Peebles potential to the DE density (cf. Eq. ( |2.17[ )) we get the 
range M ~ 5 • 10 2 — 5 eV. For reader's convenience we put the model's parameters together with 
some other quantities in Table |[ Out of curiosity we also added there an entry with an ultralight 
Dirac mass mp = 1 KeV. 

At high temperatures T > Tew when mrj = the model is in a stable phase whose parameters 
(up to a factor 2) coincide with that of the above case of a single Dirac fermion Ur\. The mass 

of active neutrino m a = 0, while the sterile neutrino is light and relativistic m s = m,R » ^(1) ■ 

/ , \2/(«+2) 

m(m/t) 

At T < Tew behavior of the model with Majorana neutrinos demonstrates some new qualitative 
features comparatively to the case [^] due to: (i) the new mass scale mrj which sets in below Tew ; 
( ii) decreasing m a with growing quintessence. The Dirac mass brings in a new critical temperature 

2 

T c d = -zm D , (3.4) 



below which the sign of the second term in the minimum equation (3.3) is always negative and no 
nontrivial solution for decreasing quintessence potentials exists. 2 

The evolution of the model below Tew towards the critical point of the first order phase tran- 
sition is qualitatively similar to that shown in Fig. [I]. The critical temperature where the maximum 
and minimum of the thermodynamic potential merge into an inflection point is found along with 
the critical value of the mass m«; 

1 + q— J ) , m c R n$ a M, (3.5) 

where the numbers ~ ~ ^"(l)- For a = 1 we find £i w 2.86 and #i w 6.42. As one can see 
from Table [j] in the physically interesting range of parameters T c T c d, so the critical temperature is 



determined in fact solely by the Dirac mass mo (compare to Eq. ( [2.1 1| )), while the quintessence at T c 



2 For large mo ~ Tew this opens an interesting opportunity for the future work to explore trial DE potentials other 
than the Ratra-Peebles one. 
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T 


T>T EW ~ 10 2 GeV 


T = T C ^ 0Am D 


T — T 


M 


m R = m s ~ m 1 y ) 


m L R « 6.42M 




m D 


m a =m D = 




m a - fixed 


5 • 10 2 eV 


m R < 10~ 3 eV 


3 KeV 


10 15 GeV 


10 2 GeV 


m s < 10~ 3 eV 


10 2 GeV 


10 15 GeV 




m a 


10 2 GeV 


10~ 2 eV 


5 eV 


m R < 10~ 6 eV 


30 eV 


10 5 GeV 


1 MeV 


m s < 10~ 6 eV 


1 MeV 


10 5 GeV 




m a 


1 MeV 


10~ 2 eV 


0.3 eV 


m R < 10~ 9 eV 


2eV 


10 8 eV 


1 KeV 


m s < 10~ 9 eV 


1 KeV 


10 8 eV 




m a 


1 KeV 


10~ 2 eV 



Table 1: Neutrino masses at different temperatures (a = 1) for three choices of the Dirac mass mo- The 
current active neutrino mass is set m a = 10~ 2 eV and the quintessence potential scale M is chosen to match 
the current DE density. All the parameters in the table are defined in the text. 



is controlled by the Ratra-Peebles scale M, similar to Eq. ( 2.13j ). The quintessence field being heavy 
at high temperatures ~ T, becomes lighter with cooling of the Universe, reaching ~ M at 
T c , and vanishes below T c . We have checked that the sound velocity and the compressibility vanish 
at T c , thus this the spinodal decomposition point and the transition occurrs via rapidly growing 
density fluctuations below T c . 



4. Discussion: Realistic Extension of the Standard Model, Dark Matter and Dark 
Energy 

It appears that the above model of Majorana neutrinos interacting with the quintessence and 
getting their masses from DE is consistent, and it encourages for further work on it. But we have 
to stress that below the phase transition the model is probed so far only by consistency checks 
between various parameters (Pde, ma, ^is) and observables. Analysis of the equation of motion 
( [2.15 ) for the quintessence field and/or by using more general methods of kinetics of the first-order 



phase transition is warranted. 

However, the model itself needs to be more realistic. It seems to us that the so-called vMSM 
[Q, ||] is the right model which has already properties needed. This is an extension of the standard 
model where three additional right-handed Majorana gauge singlets for each of the neutrino flavors 
are added. With the model's available fitting parameters it is possible [|J] to get along with the three 
light active neutrinos also the three sterile neutrinos with consistent abundances, masses, and life 
times. 3 One of the sterile neutrinos with the life time comparable to the age of the Universe has 
a mass in the KeV range and is an excellent "warm" DM candidate, while the other short-lived 
heavy neutrinos have their masses in the GeV range and are decaying DM particles. The results 
presented here allow us to anticipate that using the quintessence field for generating the right- 



handed Majorana massive terms (3.1) in this extended standard model would also allow to obtain 



In addition this model also provides a mechanism for baryon asymmetry of the Universe. 
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the observable DE density and the slow-rolling quintessence dynamics with w<p = Py/pq, « — 1. 
This study is currently work in progress. 
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